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1. 0 INTRODUCTION 



The FFTD (Fast Fourier Transform 4- Display) program can perform a Fast 

Fourier Transform or Inverse Fast Fourier Transform on 4 to 1024 real 

1 

or complex points which have been stored on a LAP6-DIAL or data LINC- 
tape or disk. The real and imaginary parts of the input or output data 
and the magnitude of the output data may be displayed on the scope via 
a moving window. Transformed data may also be stored on a DIAL or data 
LINCtape or disk. In addition, the scale of the displayed data can be 
user-modified over twelve different ranges. 

2.0 MINIMUxM HARDWARE REQUIREMENTS 
8K PDP-12B with EAE. 



3.0 OPERATING PROCEDURE 

3.1 Loading FFTD 

FFTD is a "load and go" program and is called from tape or disk by the 
DIAL command: 

-^LO FFTD, n ) 

where n is the tape (JZf-7) or disk (10-17) containing the program. A 
DIAL system tape must be on unit 0. (If a non-existent unit is addressed, 
NO is displayed on the scope. Press RETURN and issue the proper com- 
mand. ) 

At any time during program operation, FFTD may be restarted by press- 
ing the console keys; LING mode, I/O PRESET, and START 20. 

3.2 FFTD Displays 
The first display is: 

DISPLAY 1 SINGLE PRECISION FFT 

INPUT ON DIAL UNIT? Y/N 



LAP6-DIAL is hereafter referred to as DIAL. 
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Type y if the data file is on a tape or disk containing DIAL; type N 
if the file is on a data tape or disk, (A file copied from paper tape 
via PIP must be referenced as a data tape or disk.) 

The final user replies to all the scope displays are terminated by 
pressing LINE FEED. 

If the input is on a DIAL tape or disk, the second display is: 

DISPLAY 2 UNIT NUMBER^. 

FILE NAME 



Specify the unit number, JZf to 7 for tape, and 10 to 17 for disk, where 
the file is located and press RETURN. Then type the file name, which 
may be 1 to 8 characters long and must begin with a non-numeric charac- 
ter and not contain a ?, /, \ , or> . After typing the file name, press 
LINE FEED. Note that a file addressed by name on a DIAL tape or disk 
can not have a header block and must have been placed on the device 
only by the FFTD program. If a non-existent unit is requested, NO is 
displayed. To restart the program from LINCtape, press STOP, I/O PRESET, 
and START 20. The program must be reloaded from an RK8 or RFJgfS disk. 

The user is told if the file is not on the specified unit: 

DISPLAY 3 CANNOT FIND 

HIT RETURN TO CONT 

Press RETURN to bring back display 2 

If the input is on a data tape or disk, the second display is: 

DISPLAY 4 UNIT NUMBER 

BLOCK NUMBER^__ 

The unit may be any number from 0 to 7 for tape and 10 to 17 for disk. 
The block number must be an octal number from 0 to 777. If a data file 
with a header block is on a DIAL device, it may be accessed by this 
sequence (instead of the DIAL message) . The correct block number is 
the value in the DIAL index plus one. After the file has been located, 
the calculation must be specified. 
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DISPLAY 5 HOW MANY PTS^^ 

(4-1024 BY POWERS OF 2) 
REAL OR 
COMPLEX? R/C^ 

Powers of 2, from 2 to 10, are acceptable, permitting 4 to 1024 points. 
Type R if the data is real; type C if it is complex. (Refer to Sec- 
tion 4.0 for a description of data storage format.) If there is not 
enough room between the starting block number and the end of tape to 
hold the number of points specified, display 5 will reappear. 

The calculation is further specified: 

DISPLAY 6 FFT OR DISPLAY? F/D 

TRANSFORM OR 
INVERSE? T/I^ 

If the data is just to be displayed, type D and press RETURN, Then 
type T if the data has most recently been transformed or I if it has 
not been manipulated at all or has been inversely transformed. Con- 
tinue at display 7. 

The next display is : 

DISPLAY 7 OUTPUT ON DIAL UNIT? Y/N^ 

Type Y if output is to a DIAL tape or disk; type N if output is to a 
data tape or disk. 

A reply of Y to display 7 (DIAL tape or disk) causes the display: 

DISPLAY 8 UNIT NUMBER 

FILE NAME 



These answers have the same restrictions as the input display, display 
2. If there is not enough space on the DIAL tape/disk to hold the out- 
put data, the next display is: 

DISPLAY 9 NO SPACE 

HIT RETURN TO CONT 

Press RETURN to bring back display 7. 

If a file already exists with the specified name, the next display is: 
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DISPIiAY 10 REPLACE? Y/N_ 

Type Y or N to replace or not to replace the file. A reply of N will 
cause display 8 to reappear. If the file is to be replaced, but the 
new file is larger than the old file, display 9 will reappear. 

If output is to a data tape or disk, the next display is: 

DISPLAY 11 UNIT NUMBER 

BLK NUMBER^~__ 

The answers have the same restrictions as the input display, display 4. 

If there is not enough space from the starting block number to the end 
of the tape to hold the output data, display 9 will reappear* 

The program will now read in the data, perform a Fast Fourier Transform 
or Inverse Fast Fourier Transform, and write the results as complex 
data pairs onto the specified tape or disk. 

When the transform is completed or if just displays are desired, the 
following message is displayed: 

DISPLAY 12 WHICH DISPLAY? 

R(EAL) 

I (MAGINARY) 

M(AGNITUDE) 

S(CALE FACTOR) 

LINE FEED (RESTART) 

Type R, I, M, or S and LINE FEED to obtain the desired display. The 

scale factor is displayed as a decimal number (0-12) . (Refer to Sec- 
tion 6.0, Data Scaling, for an explanation of the scale factor.) (The 
magnitude, M, for a+b is M = vTa^TF^r. ) 

If the display is less than 512 points, it will be stationary and 
centered on the scope. If it contains 512 or more points, the display 
can be moved in either direction using A/D knob 0. 

A cursor which can be moved by rotating A/D knob 1 will ride along the 
curve. Associated with the cursor are four octal words displayed in 
the top left corner of the scope, one beneath the other. The first two 
words are the absolute 15-bit core address of the cursor point. The 
third word is the contents of the displayed core address, i.e., the 
actual 12-bit value in the data buffer of the data word that corresponds 
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to the cursor point* The fourth word is the scope Y coordinate of the 
cursor point. The fourth word is a relative value and depends upon the 

y scale factor and Y offset. Because the data is scaled to nine bits 
prior to display, the fourth word or Y coordinate will range from 0001 
to 1000^, where 0001 corresponds to the bottom of the scope and 1000 
to the top. 

The curve can be expanded in the Y direction by typing a 1 or decreased 
t>y typing Q. Twelve different ranges are possible. As the display is 
enlarged, no check is made against losing significant digits of large 
values because the user may wish to expand small features of the dis- 
play. Therefore, as the display is enlarged, large values may sudden- 
ly decrease in size as significant digits are lost. 

The magnitude display is shown at half scale initially. If the values 
allow, the number 1 can be typed once to show the display at full 
scale. 

Pressing RETURN will cause display 12 to reappear. As many displays 
as desired may be requested. Subsequent displays will be initially 
shown at the same range as the preceding display. Pressing LINE FEED 
without entering a character will cause display 1 to reappear. 

4.0 EXAMPLE 

This section provides examples of the displays which result from a 
transform performed on a square wave of 512 points and from an inverse 
transform performed on the resulting coefficients. 

4 . 1 Input Display 

Consider a square wave"^ of 512 real points which has the following 
format on tape or disk: 

Address Value 



The displays shown on the following pages are adaptations and are for 
demonstration purposes only. 



77 
100 




111 



1000 
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Address 


Value 




20;?0 


211 








377 


1000 




2000 


477 


1000 




0000 


511 


1000 


600 


2000 


eii 


1000 


700 


0000 


111 


1000 




If the input is displayed, there will only be a REAL display, it will 
look as follows, assuming the cursor is to the extreme left and the 
display is not moving. 



0001 
0000 
2000 
0601 



The first two values in the upper left hand corner are the address of 
the point on which the cursor is resting. When the cursor is at the 
extreme left, it indicates location 0000 of field 1. The third value 
is the contents of that memory location, in this case, 2000. The fourth 
value is the position of the cursor with respect to the bottom of the 
screen. [1 = bottom, 401 = x axis (middle), 1000 = top.] 



4.2 Transform Displays 
4.2.1, Real Display 



0001 
2000 
0000 
0401 
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Moving the cursor to the highest point in the display will change the 

value display to: 0001 

2400 
2000 
0601 



This is the DC component of the wave. 



4.2.2 Imaginary Display 



0001 
2000 

0000 
0401 



Moving the cursor to the lowest point produces the values: 



0001 
2374 
6567 
0257 



Moving the cursor to the highest point displays: 



0001 
2404 
1214 
0522 



4.2.3 Magnitude Display 



0001 
2000 

0000 
0401 



Moving the cursor to the highest point gives the following display: 

0001 
2400 
1000 
0501 



Because the magnitude of maximum values causes overflow, a factor of 2 
is removed during computation. Therefore, the values displayed are 
half scale; type the key "1" once to display the magnitude at full 
scale. 

4.2.4 Scale Factor Display 

The scale factor has a value of 1. To obtain the actual coefficients, 
rest the cursor on the desired point and shift right the third value 
of the corner display the number of bits equal to the scale factor. 
In this example, the highest value of the real display is 2000. Shift- 
ing it right by the scale factor (=1) yields 1000, the actual value of 
the DC component, which in binary is 001 000 000 000. Because the bi- 
nary point is to the right of the sign bit, the actual value is +,01^. 

4.3 Inverse Transform Displays 

The output of the transfer was 512 complex points. The inverse yields 
the following displays: 

4.3.1 Real Display 

0001 
0000 

0764 

0477 



The third value, 07 64, is a deviation from 1000, the exact value. At 
this time there are 2 scale factors involved. The relationship between 
the computed results and the original data is: 

results = [{original data) *2 factors points 

Reducing the equation for the first point yields: 

1000Q = [(2j2f0j2fg)*2^]/1000g 

2^ = 21^*2^2^ 

=2^ 
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4.3*2 Imaginary Display 



0001 
1000 
0001 
04)2(1 



The values are very small and are the result of imprecision in the com-- 
putations. 

4.3.3 Magnitude Display 

0001 
2000 
0372 
0440 



As in the magnitude display of the transform, the values displayed are 
half scale. Because the imaginary components are essentially zero, the 
magnitude, when doubled, equals the real values. 

4.3.4 Scale Factor Display 

The scale factor has a value of 7. 

5.0 DATA STORAGE 

The data must be stored sequentially on tape or disk in a binary file 
starting at the beginning of a block. If the data is complex, the 
real parts are grouped together followed by the imaginary parts, if 
any. If there are none, the program will create imaginary parts of 
value zero. The inpat and output data are in the form of binary frac- 
tions. For output data, the location following the last imaginary part 
contains the scale factor (refer to Data Scaling, Section 6.0). A 
file of complex values are stored in the following format: 
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2n 



REALi 
REAL 2 

IMAG, 
n 



SCALE FACTOR 



-only present if file is 
generated by the FFTD 
program. 



6.0 DATA SCALING 

All calculations in FFTD are done with single precision fixed point 
signed binary fractions. The binary point is located between bit 0 
and bit 1, leaving an 11 bit signed mantissa* Bit 0 is used as a 
sign bit. Negative numbers are formed by taking the two's complement 
of the positive binary fraction, so all inputs must be scaled in mag- 
nitude to les3 than one. The outputs are also formatted as above. 

In order to preserve precision, it is sometimes necessary to divide 

by 2 in a computation. As a result, a pseudo floating point format 

has been adopted in which a variable scale factor (or exponent) is im-- 

posed' on all the Fourier coefficients. This scale factor or pseudo 

exponent is found in item SCAL after each transform has been completed. 

It is also stored after the last imaginary part on tape or disk, . The 

values stored on tape or disk are the Fourier coefficients multiplied 
SCAL 

by 2 . Because in binary notation shifting a number right one bit 
is equivalent to dividing by two, to retrieve the coefficients themselves, 
shift each number right by the number of bits equal to the value of the 
scale factor. In the case of the inverse transform, the time samples 
are the values in memory multiplied by 2""^^^^. If, however, the in-- 
verse transform was performed on normalized transform data, the results 
are equal to ([ (original data) *2^]/no. of points )where n equals the sum 
of both scale factors. To retrieve the time samples, shift left each 
number by the value of the scale factor. 

7,0 SUBROUTINES USED 

Manipulation of the DIAL and data LINCtapes and disk is done using the 
program MILDRED (DEC-12-FZDA) . The question and answer displays are 
handled by QANDA (DEC-12-FISA) . The data displays are handled by DISPLAY 
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(DEC-12-FLSA) . A modification of FFTS-C (DECUS #8-144) is used to per 
form the Fourier Transforms. 



8.0 ALGORITHM DESCRIPTION 



The Fast Fourier Transformation enables computation of the power spec- 
trum of a time series in a minimum of time. Specifically, it permits 
the discrete Fourier transformation 



N-1 



X, e 
k 



•2Trijk/N 



i= >PT 



,N-1 



of a series on N equally spaced time samples (where N is a power of 2) 
The time required is proportional to N^log^N, whereas previous methods 
required times proportional to N. This gives a reduction in computa- 
tion time of l-log^N/N or over 99 percent for N=lj2(24. The algorithm 
makes use of the fact that 

to reduce the number of manipulations necessary for a transformation. 



9.0 CORE CHART 



Field 0 



SEGMENT JZf 

PAGE 0 



SEGMENT 
SEGMENT 



*4j2l0 
*1400 
1 
2 



SEGMENT 3 



Field 1 
0 

2000 



IFFT 
FFT 

DISPLAY 
MILDRED 
MONITOR 
QANDA 

Data display code 
FDV table 
RWPARM table 
Questions 
Sine Table 



Buffer - real parts 
Buffer - imaginary parts 
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10.0 PROGRAM REGION DESCRIPTION 



10.1 Routines 



IFFT 



Take the Inverse Fourier Transformation of the 
data in field 1. The results are in bit inverted 
order (refer to the SORTX routine) . 



FFT 



Take the Fourier Transformation of the data in 
field 1. The results are in bit inverted order 
(refer to the SORTX routine) • 



SORTX 



Sort the data from bit inverted order to sequen- 
tial order. Bit inversion means simply the pro- 
cess of re-ordering the bits in a binary number. 
For instance, the binary number 001 bit inverted 
is just 100 (=4). For example, to locate S^ in 
memory for a 16 point transformation (N^^'IS, n=4) , 
write 5 as a binary number of n=4 bits, 5^^^0101^, 
Then reverse the order of these bits to IJZfljZf^. 
This means S^ is stored in position 10. Physic- 
ally, then, S^ of the real parts is to be found 
in location XRTAB+9. 



MULTIP - Perform a rounded single precision signed multiply 
using EAE. The CAL+1 contains the address of the 
multiplicand. The AC contains the multiplier. 
Exit with the product in the AC. 

INVRT - Reverse the bits of the number contained in the 
AC. 



TRIGET - Fetch sine and cosine values. Specifically, if 
the AC=K on entry, the values of sin (2irK/N) and 
cos (2ifK/N) are fetched from an internal trig 
table. K must be >N/2, A register COSINE con- 
tains the cosine value and the AC contains the 
sine value on exit. 



ADDR - Perform a single precision add with rounding. 
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This subroutine generates a moving window dis- 
play with a cursor riding on the curve. For 
more information refer to the DISPLAY document, 
DEC~12-FLSA-D. 

Display the question: FROM DIAL UNIT? Y/N_ 
If the answer is Y, jump to UNTFIL; if N, jump 
to DATTAP; if neither, redisplay the question. 

Jump to the subroutine ASK2 to display: 

UNIT NUMBER 
FILE NAJ^dE 



If the unit number is illegal, jump to ASK2 
again to redisplay the question. If legal, jump 

to LOOKUP with the address of the File Descrip- 
tion Vector (hereafter referred to as FDV) para- 
meter list in the AC. If the file cannot be found, 
display the message: 

CANNOT FIND 

HIT RETURN TO CONT 

When RETURN is hit, jump back to UNTFIL. If the 
file is found^ jump to MOVINP. 

Jump to the subroutine ASK3 to display: 

UNIT NUMBER 
BLK NUMBER^^r 

If an illegal value is entered, jump back to 
DATTAP. If all the input is legal, fall through 
to MOVINP. 

Jump to FDV2RW to move the input information from 
the FDV to the read/write parameter list. Fall 
through to PTS. 

Display: NUMBER OF PTS 

(4-1024 BY POWERS OF 2) 
REAL OR 
COMPLEX? R/C 
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Set Bl to the address of the answer buffer, 
MPLIER to 12 and UPLEGL to -71 (-9) because 
the number of points is entered as a decimal 
value. Set the AC to the largest legal value, 
2000, and jump to CONV. If the answer is an 
illegal value jump back to PTS; store the value 
in N and store its I's complement in TEMPI. 
Since the number of points must be an integral 
power of 2, only one bit in TEMPI may be set. 
Bit 11 is the exception to one bit being a power 
of 2. Check bit 11 first, then rotate the value 
adding up the number of bits set. If the total 
is not 1, jump back to PTS. Otherwise fall through 
to ROTl. 

Compute the power of 2 by rotating right the 
value in TEMPI and stepping B2 until the bit 
that is set is encountered in bit 11, Fall 
through to STAMU, 

Store the power of 2 in NU, If the power is 
less than 2, jump back to PTS. Otherwise load 
the AC with the number of points *2 and jump to 
NUMBKS to compute the number of blocks needed 
to hold the output. Store the value in FDV+7. 
Store it also in RWPARM+3 since, for complex 
data, the input and output data consist of the 
same number of blocks. If the answer to the 
second question is not R; jump to IFCOM. If it 
is R, the input consists of half as many words 
as the output. Load the AC with the value of N 
and jump to NUMBKS to compute the number of in- 
put blocks. Store the value in RWPARM+3. Set 
REALPG and jump to CKEND. 

If the answer is C, clear REALFG and fall through 
to CKEND. Otherwise jump back to PTS to redisplay 
the question. 

.If there is not enough room between the starting 
block number and the end of tape to hold the num- 
ber of points specified, jump back to PTS. If 
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the number of output words is 400 or greater, 
another block will be needed to hold the scale 
factor • Increment FDV+7. Fall through to IFFFT. 



IFFFT 



Display : 



EFT OR DISPLAY? 
TRANSFORM OR 
INVERSE? T/I 



F/D 



If the answer to the first question is D, set 
DISFLG to indicate that the data will only be 
displayed. If F, clear DISFLG to indicate that 
a Transform or Inverse Transform will be performed. 

If the answer to the second question is T, clear 
FTFLG; if I, set it. If DISFLG is set, jump to 
DISPLY to display the data. Otherwise, jump to 
OUTQES , 



OUTQES 



Display the question: OUTPUT ON DIAL UNIT? Y/N_ 
If the answer is Y jump to OUTUNT; if N jump to 
ONDAT; otherwise redisplay the question, 



OUTUNT 



Jump to the subroutine Aa<2 to display; 



UNIT NUMBER 
FILE NAME 



If an illegal value is input, redisplay the ques- 
tion. Otherwise jump to ENTER with the address of 
the parameter list in the AC. If a file with the 
specified name already exists, jump to SAMNAM. If 

there is not enough space to hold the output data, 
jump to NOSPAC. If it is a new file and there is 
enough space "to hold it, fall through to RDDATA. 

RDDATA - Clear 4000 words of field 1 and read in the input 
data. If REALFG is 0, the data is complex - move 
the imaginary parts to start at location 2000. If 

it is non-zero, the data is real and nothing need 
be done. Jump to PROC. 



PROC 



If IFTFLG is 0, jump to FT to do a Transform. 
Otherwise, fall through to do an Inverse Transform. 
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Jump to the subroutine IFFT to do an Inverse 

Transform on the input data. Then jump to the 
subroutine SORTX to sort the coefficients into 
sequential order from bit inverted order. Jump 
to STSCAL to store the scale factor which is equal 
to NU-SCAL, The data should be shifted by this 
value. 

Jump to the subroutine FFT to transform the input 
data. Then jump to the subroutine SORTX to sort 
the coefficients into sequential order from bit 
inverted order. The scale factor is the value in 
SCAL and equals the number of bits by which the 
data should be shifted right. Fall through to 
STSCAL. 

Store the scale factor in the word following the 
last imaginary part. Move the imaginary parts from 
200 to immediately behind the real parts. 



Jump to the subroutine FDV2RW to move the output 
parameters from the FDV to the read/write parameter 
list. Write the data onto the output tape and jump 
to DISPLY. 

Jump to the subroutine ASK to display the message: 
NO SPACE 

HIT RETURN TO CONT 

When RETURN is hit, jump to OUTQES . 

Jump to the subroutine ASK to display: 

ALREADY EXISTS 
REPLACE? Y/N_, 

If the answer is Y, jump to REPL; if it ia N, jump 
to OUTUNT. If it is neither, redisplay the ques- 
tion. 
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REPL 



Try to replace the existing file with the new 
file. If the new file is longer, jump to NOSPAC, 
If the replacement is successful, jump to RDDATA. 



ONDAT 



Jump to the subroutine ASK3 to display: 



UNIT NUMBER 
BLK NUMBER " 



If an illegal value is entered^ redisplay the 
question. If there is not enough space between 
the specified block number and the end of tape to 
hold the output data, jump to NOSPAC. Otherwise, 
jump to RDDATA. 



10.2 Subroutines 



FDV2RW 



Transfer the unit number, starting block niamber^ 
and number of blocks from the FDV parameter list 
to the READ /WRITE parameter list. 



NUMBKS - Enter with the number of words in the AC, Con- 
vert this value to blocks by counting the number 
of times 400 can be subtracted from it before the 
value becomes negative • Return with the number 
of blocks in the AC. 



ASK2 



Jump to OCTL to set MPLIER to 10 and UPLEGL to 
-67 (-7) because the unit number is input as an 
octal number. 



Display: UNIT NUMBER 
PILE NAME 



by jumping to the subroutine ASK with the ad- 
dress of QUES2 in the AC. Set Bl to the address 
of the answer buffer and jump to the subroutine 
CONV with the largest legal unit number, 17, in 
the AC. If the value is illegal, return to CALL+1, 
If legal , store it and the file name in the FDV 
parameter list. Fill the file name out to 8 
characters with 77»s. Return to CALL+2. 
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Display: UNIT NUMBER 
BLK NUMBER 



by jumping to the subroutine ASK with the address 
of QUES3 in the AC. Set Bl to the address of the 
answer buffer and jump to OCTL to set MPLIER to 10 
and UPLEGL to -67(7) because the unit and block 
nximbers are input in octal. Jump to subroutine 
CONV with the largest legal unit number^ 17, in 
the AC. If the value is illegal, return to CALL+1. 
Otherwise, store it in word 0 of the FDV parameter 
list. Bl is now pointing to the block number. 
Jump to CONV with the largest legal block number, 
777, in the AC, If the value is illegal/ return 
to CALL+1, If legal, store it in word 6 of the 
FDV parameter list. Return to CALL+2. 

CONV is entered with the largest legal value in 
the AC and Bl pointing to the address - (1 half 
word) of the first character to be converted. 

Store the I's complement of the largest legal value 
in TEMP2 and clear TEMPI. UPLEGL contains a -71 (-9) 
or -67 (-7) and MPLIER contains a 10 or 12 depending 
on whether the number to be converted is in decimal 
or octal. Extract a character and compare it against 
an ASCII 0 and the contents of UPLEGL, If it is a 
legal value, jump to MULPLY which will multiply the 
value in TEMPI by the contents of MPLIER and add 
the digit being converted to it. Repeat the pro- 
cedure until a character is found which is not be- 
tween 0 and UPLEGL. If it is not a 34, 74, or i3, 
it is an illegal character: return to CALL+l. A 
34 or 74 indicates the end of the input field; a 
0 indicates the end of the input. Compare the con- 
verted value in TEMPI against the maximum legal value 
in TEMP2. If the value is legal return to CALL+2; 
otherwise return to CALL+1. 

OCTL sets MPLIER to 10 and UPLEGL to -67 (-7) so 
that CONV will convert an octal number. 
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ASK is entered with the address of the display in 
the AC. Store it in the parameter list and jump 
to QAINIT to display the message. Refresh the 
display until the answer is input. Return to the 
calling routine. 

This region is entered either after the Transform 
or Inverse Transform is completed or in response 
to a D in answer to the display: FFT OR DISPLAY? 
F/D_. Since the data is manipulated in prepara- 
tion for each display it must be read in before 
each display. After reading in the data, display: 

WHICH DISPLAY? 

R(EAL) 

I (MAGINARY) 

M(AGNITUDE) 

S(CALE FACTOR) 

LINE FEED (RESTART) 

If the answer buffer contained 0, just LINE FEED 
was hit: jump to IFDIAL to restart the program. 
Otherwise jump to WCHDIS. 

Jump to DPIMAG, DPMAG, DP REAL, or DPSCAL if the answer 
was I, M, R, or S, respectively . Otherwise redisplay 
the question. 

If REALFG is non-'zero, the input is real and no 
Transform was performed. Therefore, there are no 
imaginary parts to display; redisplay the question. 
If REALFG is zero, check IPTFLG. If it equals zero, 
either an Inverse Transform was performed or the 
original data is just being displayed. In either 
case the data is in the right order. If IFTFLG 
is non-zero, a transform was performed. The posi- 
tive half of the curve is first followed by the 
negative half and the signs are reversed. Swap 
the halves and reverse signs before jumping to 
PREPAR. 

Check IFTFLG for the same reason as in DPIMAG. 
The only difference is that the signs of the real 
parts are not reversed. 



19 



If less than 1000 points are to be displayed, 
the display will not move and the points dis- 
played will be centered on the scope. To achieve 
this, LEFTX is set to the l*s complement of 
-10jaf0+(lj2fj30=# of points) /2,MINPTS to the 2's 
complement of the number of points, and MVDIS 
to the instruction CLR. Jump to SHOWIT. 

If 1000 or more points are to be displayed, the 
display will fill the scope and will move. To 
achieve this, LEFTX is set to the I's complement 
of 1000, MINPTS to the 2's complement of 1000 
and MVDIS to the instruction SCR 4. Fall through 
to SHOWIT. 

Jump to the subroutine IDORA to display the data. 
The six parameters following the call to IDORA are 
in order: the memory field of the lower address, 
the lower address, the memory field of the higher 
address, the higher address, the Y offset of the 
display and the scale factor of the data. Both 
fields are always 1, the lower address is always 
0. The higher address is set in the region DISPLY. 
The Y offset is always 0; therefore the baseline is 
half way up the scope. The scale factor is the 
instruction SCR plus the nxuhber of bits to scale 
the data right before displaying it. Since IDORA 
displays only the right nine bits, if the left 
three bits are significant, the data must be scaled 
right three before displaying it. 

Jump to RDORA to refresh repeatedly the display 
until a key on the teletype is hit. If the RETURN 
is hit, jump to REDPLY which jumps to DISPLY to 
redisplay the question: WHICH DISPLAY? If a 1 
is entered, jump to LARGER to blow up the display. 
If a Q is hit, jump to SMALLR to decrease its size. 
If anything else is entered, ignore it. 

If the instruction at SIZE contains a shift of 11 
bits, a bigger shift would be meaningless. Jump 
back to RFRSH. Otherwise, increment the value of 
the shift and jump to SHOWIT. 
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LARGER - If the instruction at SIZE contains a shift of 0 
bits/ jump back to RFRSH. Otherwise decrement the 
value of the shift and jump to SHOWIT. 

DPSCAL - If REALFG is non-zero, only real parts are present, 
meaning this program did not create the file and 
therefore *there is no scale factor. Return to 
DISPLY to redisplay the question. If REALFG is 
0, the scale factor is stored after the last imaginary 
part. Convert it to ASCII decimal and display it. 

If REALFG is non-zero, the input data is real and 
no transform was performed; therefore the magnitude 
is the same as the real points. Redisplay the ques- 
tion: WHICH DISPLAY? Otherwise move the imaginary 
parts to location 2000. Set RELPTR and IMGPTR, which 
contain the effective address of the multipliers, to 
6000 since the data begins at location 0 of their 
respective segments and is fractional. Fall through 
to NXTMAG. 

NXTMAG - Square a real part and store it. Square the imaginary 
part, add the square of the real part to it, jump to 

the subroutine SQRT to get the square root of the s\im 
and store it in place of the real part. Repeat the 
process for each point. Then jump to SHOWIT to dis- 
play the magnitude. 

MDVPTS - The subroutine MOVPTS moves values from one buffer 
(address -1 in 10) in field 1 to another (address 
-1 in 11). If CMPFLG equals 1, the values are com- 
plemented as they are moved. TEMPR contains the 
2's complement of the number of values to move. 

The subroutine MVRLMG is lised to swap the first 
and second halves of the real or magnitude values. 

In the process they are moved from the buffer start- 
ing at location 0 to the one starting at 2000. 

The File Descriptor Vector parameter list is used 
by the LOOKUP, ENTER, and REPLACE sections of MILDRED. 
Word 0 contains the unit number, words 1-4 contain 



DPMAG 



MVRLMG - 
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the file name^ word 5 contains a 2 indicating 
the file is binary, word 6 is the starting block 
nuinber, and word 7 is the number of blocks. Word 
6 is filled by LOOKUP, ENTER and REPLACE. Word 7 
is filled by LOOKUP but must be supplied for ENTER 
and REPLACE. 

The Read/Write parameter list is used by the READ 
and WRITE sections of MILDRED. Bits 0-2 of word 0 
contain the field, bits 9-11 contain the unit. 
Word 0 contains the starting address, word 1 the 
starting tape block number and word 2 the number 
of blocks. 

The subroutine SQRT is entered with a value in the 
double precision location DPSQ. It returns with 
the square root in the AC. 
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10. 3 Symbols 



N Nuitiber of words in computation 

NU Power of 2 of value of N 

L Index to show what array is being constructed 

S Gives spacing between node pairs in the Lth array 

N0VER4 Storage for N/4 

MAXNU Power of 2 of largest table size (13) 

MN0VR2 Storage for N/2 

QR Pointer to real part of X(Q) 

QI Pointer to imaginary part of X(Q) 

PR Pointer to real part of X{P) 

PI Pointer to imaginary part of X(P) 

Q Numerical index Q (=0,1, . . . ,N-1) 

P Numerical index P (=0, . , • • • #N-"1) 

K Number in the node being operated on 

C Interrupts computation of Lth array every S passes- 

ADD2 Used by subroutine ADDR as data (addend) 

Used by monitor as a temporary location 
TEMPR Temporary storage register for real parts 

Used by monitor as a temporary location 
SINE Temporary storage for sin (S*PI*K/N) 

Used by monitor as a temporary location 
COSINE Temporary storage for cos (2*PI*K/N) 

Used by monitor as a temporary location 
GR Real part of product (W *X(P)) - temporary storage 

Used by monitor as a temporary location 
GI Imaginary part of product (W *X(P)) - temporary storage 

SCAL Pseudo exponent of Fourier coefficients 

SHFLAG If =1, add with shift; if =0, add without shift 

SHFCHK Indicates if all X*s in an iteration are <.5 

DISFLG If 7^0, the data will just be displayed 

IFTFLG If ^0, an Inverse Transform was performed 

REALFG If ^0 , the data does not contain imaginary parts 

DPSQ Used to save the double precision squares of the real 

and imaginary parts during calculation of the magni- 
tude. 

CMPFLG If =1, the subroutine MOVPTS will complement the 



values as it moves them 
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10. 4 Beta Registers 



Beta registers 1, 2, and 3 are used by the monitor in ASK2 and ASK3 as 
temporary pointers and counters. QANDA and MILDRED make more extensiv 
use of the Beta registers. 



11.0 ASSEMBLY INSTRUCTIONS 



The FFTD program is assembled in three sections by assembling and sav- 
ing each, then adding them together . The entire command sequence is: 



->AS MILQAN,n J 
•^SE MILQAN,nJ 
-^AS SIN256,n^ 
->SB SIN256,nJ 
">AS FFTC-lJ 
•^SB FFTC-iy 
->ZEj 

->AB MILQAN,n^ 
^AB SIN256,n J 
->AB FFTC-l,n3 
->SB FFTD,n,Lj 



where n is the unit 
containing the program 

(FFTC-1 chains to FFTC-2) 



(saves the whole program) 



12.0 SYSTEM FLOWCHARTS 



(Attached) 



13.0 PROGRAM LISTING 



(Attached) 
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^ START 20 ^ 



READ IN DIAL I/O 
ROUTINES 



DISPLAY: 



SINGLE PRECISION FFT 
, INPUT ON DIAL TAPE? Y/N 




''DISPLAY: 

UNIT NUMBER 
FILE NAME 




DISPLAY: 
UNIT NUMBER 
BLK NUMBER 




^I40VINP ^ 



MOVE INPUT PARAMETERS 
TO R/W PARAMETER TABLE 



A^Ts) 



DISPLAY: 
CANNOT FIND 
HIT RETURN TO CONT/ 




DISPLAY: 

HOW MANY PTS? 

(4-102 4 BY POWERS OF 2) 
REAL OR 
COMPLEX? R/C 



T 
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COMPUTE NUMBER 
OF OUTPUT BLOCKS 





4 



ADD 1 BLOCK 

FOR SCALE FACTOR 



4- 



N 



DISPLAY: 
FFT OR DISPLAY? F/D 
TRANSFORM OR 
INVERSE? T/I 




IFTFLG=j2f 



^ DISFLG=NON-ZERO 
I 



IFTFLG 

= NON-ZERO 



^ IFDISP 




DISPLAY: 
OUTPUT ON 
DIAL UNIT? Y/N 



^MOVE OUTPUT PARAMETERS 
TO R/W PARAMETER TABLE 
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RDDATA 



CLEAR DATA 
BUFFER 



READ IN 
DATA 




SCALE FACTOR= 
SCAL 



WORE 

P0I^ 


) AFTER LAST 
fT=SCALE FACTOR 


\ 


i 



^ NOWSTR ^ 



NOWSTR^ 



PACK IMAG PARTS 
BEHIND REAL 



JikL 



MOVE OUTPUT PARAMETERS 
TO R/W PARAMETER TABLE 



WRITE OUT 
DATA 



iDISPLY 
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READ IN 
DATA 




DISPLAY: 
WHICH DISPLAY? 
R(EAL) 
I (MAGI NARY) 
M(AGNITUDE) 
S(CALE FACTOR) 
LINE FEED (RESTART) 








CONVERT £ 
FACTOR TC 
DECIMAL 


>CALE 
) ASCII 


s 


/ 
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UST 
AL PART? 
7 



YES 



DISPLl 




NO 



MOVE IMAG PARTS 
TO 2j3i2fj2f 



SQUARE A 
REAL PART 




SQUAR 
IMAG 


E AN 
PART 






ADD SQUARES 


\ 





^SQUARE ROOT ^ 





/ 


STORE 
OF REi 


IN PLACE 
\L PART 
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YES 



LEFTX = "IW+djZfjajaf-* OF PTS)/2 (I'S COMP) 
MVDIS = CLR 
MINPTS = COSINE 
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INDEX 



Algorithm, 11 
A/D Knob 0, 4 
A/D Knob 1, 4 



Block Number, 2, 4, 14, 17, 18 



Core Chart, 11 
Cursor, 4, 5, 6 



Data Display, 3, 4, 5, 6, 7, 8, 

9, 19, 20 
Data Scaling, 5, 10 
Data Storage, 9 

Disk, 1, 2, 3 
DISPLAY, 10, 13 



Examples, 5 

Imaginary Displays, 7, 9 
Input Display, 5 
Magnitude Displays, 7,9 
Real Displays, 6, 8 
Scale Factors, 8, 9 
Transform Display, 6 



File Format, 9 

File Name, 2, 3, 13, 15 



Hardware Requirements, 1 
Header Block, 2 



Input Data, 1, 2, 3, 9, 12, 14, 

la. . 

inverse Transform, 1, 3, 8, 15 



LINCtape, 1, 2, 3 
Loading Procedure, 1 



Magnitude, 1, 4, 5, 19, 21 
MILDRED, 10 

Moving Window, 1, 4, 13, 20 



Octal Words, 4 

Output Data, 1, 3, 9, 14, 15 



Points 

Complex, 1, 3, 4, 8, 9, 14, 
15 



Points (cont'd) 

Imaginary, 1, 4, 9, 15, 16, 19, 
21 

Real, 1, 3, 4, 6, 9, 14, 15, 19, 
21 

Program Routines, 12 
Program Subroutines, 17 
Program Symbols, 23 



QANDA, 10 



Replace File, 4, 16 
Restart Procedure, 1, 2 



Scale Factor, 4, 5, 8, 9, 16, 19, 
20 



Terminators ,2,5 
Transform, 1, 3, 6, 11, 15 



Unit Number, 2, 3, 4, 13, 15, 17, 
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READER^S COMMENTS 



FFTD 



Digital EquFpmenf Corporation maintains a continuous effort to improve the quality and usefulness 
of Its publications. To do this effectively we need user feedback — your critical evaluation of 
this manual . 

Please comment on this manual's completeness, accuracy, organization, usability, and read- 
ability. 



Did you find errors in this manual? If so, specify by page . 



How can this manual be improved? 



Other comments? 



Please state your position . Date: 

Name: _ _ Organization: 



S^J'eet: — Department: 

— State:- Zip or Country. 



INDEX 



Algorithm, 11 
A/D Knob 0, 4 
A/D Knob 1, 4 



Block Number, 2, A, 14, 17, 18 



Core Chart, 11 
Cursor, 4, 5, 6 



Data Display, 3, 4, 5, 6, 7, 8, 

9, 19, 20 
Data Scaling, 5, 10 
Data Storage, 9 
Disk, 1, 2, 3 
DISPLAY, 10, 13 



Examples, 5 

Imaginary Displays, 7, 9 
Input Display, 5 
Magnitude Displays, 7,9 
Real Displays, 6, 8 
Scale Factors, 8, 9 
Transform Display, 6 



File Format, 9 

File Name, 2, 3, 13, 15 



Hardware Requirements , 1 
Header Block, 2 



Input Data, 1, 2, 3, 9, 12, 14, 

19 

Inverse Transform, 1, 3, 8, 15 



LINCtape, 1, 2, 3 
Loading Procedure, 1 



Magnitude, 1, 4, 5, 19, 21 

MILDRED, 10 

Moving Window, 1, 4, 13, 20 



Octal Words , 4 

Output Data, 1, 3, 9, 14, 15 



Points 

Complex, 1, 3, 4, 8, 9, 14, 

15 



Points (cont'd) 

Imaginary, 1, 4, 9, 15, 16, 19, 
21 

Real, 1, 3, 4, 6, 9, 14, 15, 19, 

21 

Program Routines, 12 
Program Subroutines , 17 
Program Symbols, 23 



QANDA, 10 



Replace File, 4, 16 
Restart Procedure, 1, 2 



Scale Factor, 4, 5, 8, 9, 16, 19, 
20 



Terminators ,2,5 
Transform, 1, 3, 6, 11, 15 



Unit Number, 2, 3, 4, 13, 15, 17, 
18 



HOW TO OBTAIN SOFTWARE INFORMATION 



Announcements for new and revised software^ as well as programming 
notes, software problems, and documentation corrections are published 
by Software Information Service in the following newsletters. 

Digital Software News for the PDP-8 & PDP-12 

Digital Software News for the PDP-ll 

Digital Software News for the PDP-9/15 Family 

These newsletters contain information applicable to software available 
from Digital's Program Library, Articles In Digital Software News up- 
date the cumulative Software Performance Summary which is contained 
in each basic kit of system software for new computers. To assure that 
the monthly Digital Software News Is sent to the appropriate software 
contact at your installation, please check with the Software Special- 
ist or Sales Engineer at your nearest Digital office. 

Questions or problems concerning DEC software should be reported to 
the Software Specialist. In cases where no Software Specialist Is avail- 
able, please send a Software Performance Report form with details of 
the problem to; 

Software Information Service 
Digital Equipment Corporation 
146 Main Street, Bidg. 3-5 
Maynard, Massachusetts 01754 

These forms which are available without charge from the Program 
Library, should be fully filled out and accompanied by Teletype output 
as well as listings or tapes of the user program to facilitate a complete 
investigation. An answer will be sent to the individual and appropriate 
topics of general interest will be printed in the newsletter. 

New and revised software and manuals. Software Performance Report 
forms, and software price lists are available from the Program Library. 
When ordering, include the document number and a brief description of 
the program or manual requested. Revisions of programs and documents 
will be announced in the newsletters. Direct all inquiries and requests 
to: 

Program Li brary 
Digital Equipment Corporation 
146 Main Street, Bidg . 1-2 
Maynard, Massachusetts 01754 

Digital Equipment Computer Users Soc ie t y ( D ECU S) maintains a user 
library and publishes a catalog of programs as well as the DECUSCOPE 
magazine for its members and non-members who request it. For further 
information please write to: 

DECUS 

Digital Equipment Corporation 
146 Main Street, Bidg . 3-5 
Maynard, Massachusetts 0 1 754 
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